Vibration State Identification of Hydraulic Units Based on Improved Artificial Rabbits Optimization Algorithm

To improve the identification accuracy of the vibration states of hydraulic units, an improved artificial rabbits optimization algorithm (IARO) adopting an adaptive weight adjustment strategy is developed for optimizing the support vector machine (SVM) to obtain an identification model, and the vibration signals with different states are classified and identified. The variational mode decomposition (VMD) method is used to decompose the vibration signals, and the multi-dimensional time-domain feature vectors of the signals are extracted. The IARO algorithm is used to optimize the parameters of the SVM multi-classifier. The multi-dimensional time-domain feature vectors are input into the IARO-SVM model to realize the classification and identification of vibration signal states, and the results are compared with those of the ARO-SVM model, ASO-SVM model, PSO-SVM model and WOA-SVM model. The comparative results show that the average identification accuracy of the IARO-SVM model is higher at 97.78% than its competitors, which is 3.34% higher than the closest ARO-SVM model. Therefore, the IARO-SVM model has higher identification accuracy and better stability, and can accurately identify the vibration states of hydraulic units. The research can provide a theoretical basis for the vibration identification of hydraulic units.


Introduction
With the continuous development of the new energy industry, hydropower, a type of renewable and clean energy, has received more and more attention. As the core equipment of hydroelectric power generation, the operating states of hydraulic units are related to power generation efficiency, economic benefits and production safety [1]. Vibration is an important factor that affects the stable operation of hydraulic units. Comprehensive, effective and accurate classification and identification of vibration states are the premise of monitoring the operating states of hydraulic units, and it is very necessary to improve the identification accuracy of the vibration states of hydraulic units to judge the stable operation of hydraulic units.
Since the hydraulic unit is subjected to the combined action of mechanical, hydraulic and electrical factors, its vibration signal shows obvious non-stationary and non-linear characteristics [2]. For the identification of the vibration states of hydraulic units, processing and extracting the feature vector of the vibration state signals is an indispensable effective information. Therefore, the VMD method is used to decompose the different states of the vibration signals of the hydraulic units, select the IMF component, and extract the time-domain feature indicator of the IMF component. Through multi-source information fusion, multiple time-domain feature vectors are combined into a new multi-dimensional feature vector and finally input into the IARO-SVM model for vibration state identification. The IARO-SVM model is compared with the ARO-SVM model, the ASO-SVM model, the PSO-SVM model and the WOA-SVM model. The results show that the IARO-SVM model has a higher identification accuracy and can effectively distinguish different signal states of vibration states.
The main motivation of this paper is summarized as follows.
• The stable operation of hydraulic units is very important for the safe production of hydropower stations. Vibration state monitoring is the premise of analyzing whether the hydraulic unit is operating stably. It is very necessary to improve the identification accuracy of the vibration states of hydraulic units to guarantee the secure operation of the hydraulic units; • Optimization algorithms are widely used in signal classification and identification. However, many optimization algorithms have some shortcomings. Therefore, according to the shortcomings of the ARO algorithm, this paper proposes an improvement strategy based on the adaptive adjustment of the weight to improve the identification accuracy of the vibration states of the hydraulic units; The main contributions of this study are highlighted as follows.
• A method of adaptively changing the inertia weight according to the current rabbit population distribution is proposed, which is used to improve the optimization ability of the ARO algorithm; • A total of 23 benchmark functions are used to test the performance of the IARO algorithm. The experimental results show that the IARO algorithm has good exploration and exploitation abilities; • Through the application of vibration state signals of the hydraulic units, the identification method of SVM optimized by the IARO algorithm is verified. The experimental results show that the IARO-SVM model can well identify the different vibration states of hydraulic units, and provide a guarantee for the stable operation of hydraulic units.
The residual section of this paper is as follows. Section 2 describes the basic principles of the VMD algorithm, the ARO algorithm and SVM. In Section 3, a method of adaptively changing the inertia weight according to the current rabbit population distribution is proposed to improve the original version. Section 4 describes testing the performance of the IARO algorithm using 23 benchmark functions. Section 5 verifies the practicality of the IARO-SVM model in engineering practice by identifying the vibration states of hydraulic units. Section 6 introduces the conclusion and future development direction.

Variational Mode Decomposition Algorithm
To solve the problems of modal aliasing and endpoint effects in signal decomposition by classical recursive decomposition algorithms such as local mean decomposition (LMD) and EMD, Dragomiretskiy et al. [27] proposed a completely non-recursive VMD algorithm. The VMD [28,29] has strong noise robustness, by which non-stationary sequences with high complexity and strong nonlinearity can be decomposed into several relatively stationary subsequences with different frequency scales [30]. Its core is the construction and solution of variational problems, which describe the VMD variational problem as follows: assume that the original signal y(t) consists of K eigenmode functions {u k (t)}, (k = 1, 2, · · · K) with different center frequencies and limited bandwidth, and constrain the sum K ∑ k=1 u k (t) of each eigenmode function to restore the original signal y(t). The constrained variational problem is expressed as follows: where δ(t) is the pulse function, j represents imaginary number; ⊗ represents convolution, y(t) is the original signal, ω k is the center frequency of each modal component u k (t), and ∂ t represents the derivation of time.
To solve the above variational model Equation (1), the quadratic penalty factor and the Lagrange multiplier are introduced to transform the above constrained variational problem into an unconstrained variational problem: where · represents the inner product, is the second penalty, α reduce the influence of Gaussian noise on signal reconstruction, and λ(t) ensure the strictness of the constraints. On the basis of the transformation variational problem, the alternating direction method of multipliers (ADMM) is used to alternately update u n+1 k , ω n+1 k , and λ n+1 to solve the 'saddle point' of Equation (2), that is, to solve the optimal solution of the variational problem. The iteration process is stopped when the following conditions are met: where ε is the discrimination accuracy. After the end of the cycle, K modal components u k (t) can be obtained by inverse Fourier transform of u k (ω).

Artificial Rabbits Optimization Algorithm
The ARO is a new bionic heuristic algorithm [31] inspired by the famous Chinese proverbs "Rabbits don't eat the grass on the edge of the nest" and "A cunning rabbit has three burrows". Figure 1 depicts the three burrows of a cunning rabbit. The ARO algorithm adopts the detour foraging and hiding strategies. Assume that each rabbit in the population has its own area, each area has edible grass and d burrows, and rabbits always visit the foraging position randomly. The mathematical model of the detour foraging is as follows [31]: Illustrates that a cunning rabbit has three burrows [31]. Assume that each rabbit in the population has its own area, each area has edible grass and d burrows, and rabbits always visit the foraging position randomly. The mathematical model of the detour foraging is as follows [31]: where v i (t + 1) is the candidate position of the ith rabbit at the time t + 1, x i (t) is the position of the ith rabbit at the time t, n is the number of rabbits, d is the dimension, T is the maximum number of iterations, • is the ceiling function. L is the step size. n 1 is subject to the standard normal distribution. c is the mapping vector.

Random Hiding (Exploitation)
To reduce the probability of being hunted, a rabbit randomly chooses a burrow to hide. The ith rabbit is in the ith burrow, and the mathematical expression is as follows [31]: To simulate this random hiding behavior of rabbits, the following mathematical equations are given [31]: where b i,r is a randomly selected burrow from d burrows. r 4 and r 5 are two random numbers between (0, 1). According to Equation (14), the ith rabbit randomly selects a burrow from its d burrows to update its position.
After the detour foraging and random hiding are completed, the position of the ith rabbit is updated as follows [31]: An energy factor needs to be designed to simulate the transition from exploration to exploitation throughout the iteration. The energy factor is defined as follows [31]: The behavior of A over 1000 iterations is depicted in Figure 2. In Figure 2, when the energy factor A(t) > 1, rabbits tend to randomly explore the areas of different rabbits for foraging in the exploration phase; thus, detour foraging occurs. When the energy factor A(t) ≤ 1, rabbits tend to randomly exploit their own burrows in the exploitation phase; thus, random hiding occurs.
where , i r b is a randomly selected burrow from d burrows. 4 r and 5 r are two random numbers between ) 1 , 0 ( . According to Equation (14), the i th rabbit randomly selects a burrow from its d burrows to update its position.
After the detour foraging and random hiding are completed, the position of the i th rabbit is updated as follows [31]:

Energy Shrink (Switch from Exploration to Exploitation)
An energy factor needs to be designed to simulate the transition from exploration to exploitation throughout the iteration. The energy factor is defined as follows [31]: The behavior of A over 1000 iterations is depicted in Figure 2. In Figure 2, when the energy factor ( ) 1

> t A
, rabbits tend to randomly explore the areas of different rabbits for foraging in the exploration phase; thus, detour foraging occurs. When the energy factor ( ) 1

≤ t A
, rabbits tend to randomly exploit their own burrows in the exploitation phase; thus, random hiding occurs.

Multi-Classification Design of SVM
SVM is a class of generalized linear classifiers for binary classification of data according to supervised learning, and its decision boundary is the maximum margin hyperplane for solving learning samples [32,33].
where ω is the weight. C is the penalty factor. ξ is the relaxation factor. b is the bias constant.
The optimal classification decision function obtained is When the order of the polynomial kernel is not 1, a SVM can be obtained. The SVM optimization problem is The optimal classification decision function obtained is The binary classification SVM can be extended to multi-classification identification based on a multi-classification strategy. The "one-to-one" (OVO) strategy is combined with SVM to form a multi-classifier, which can obtain better classification performance [34]. When solving the multi-classification problem of K-class signals, constructing a two-class SVM sub-classifier between any two categories of samples, a total of K(K − 1)/2 subclassifiers need to be constructed. In this study, since the Gaussian radial basis function (RBF) has fewer parameters, it is used as the kernel function of the SVM [17].

Improvement of Artificial Rabbit Optimization Algorithm
The ARO has the advantages of simple operation, few parameters to be set and strong optimization performance, but there are also points that can be improved, such as convergence accuracy, convergence speed, and local optima avoidance. Using dynamic inertial weights, in the early stage of algorithm iteration, due to the relatively dispersed population, larger weight values can be assigned to accelerate the global search ability of the algorithm. In the later stage of iteration, the algorithm can adaptively change the size of the weight value according to the distribution of individuals in the current population, so that it can search finely around the optimal solution and accelerate the convergence speed.
According to the inspiration of Kong et al. [35] on the improvement of the whale optimization algorithm (WOA), this paper proposes a method to improve the ARO algorithm, namely, the ARO with adaptive weight (IARO).

Adaptive Inertia Weight
The inertia weight is an important parameter in improving the ARO algorithm. The appropriate weight value can improve the optimization ability of the algorithm. However, the improper selection of linear inertia weight adjustment strategies will affect the convergence speed of the algorithm. Therefore, this paper proposes a method to adaptively change the weight value according to the current rabbit population distribution, which is as follows: where t represents the number of iterations of the current population. P iworst and P ibest are the position vector of the worst rabbit and the position vector of the optimal rabbit in the current rabbit population, respectively. The obvious difference between the proposed inertia weight and the weight in [35] is that our method uses the current population information to update the lower and upper boundaries to adaptively adjust the search space; another difference between them is that our method uses two random coefficients to compromise the maximum individual distance and the distance in each dimension. The adaptive adjustment for the weight of the current rabbit in the random hidden update position is as follows: After introducing the adaptive adjustment weight strategy, the algorithm can adaptively change the weight size according to the distribution of the current rabbit population. In the early stage of algorithm iteration, if the rabbit population falls into the local optimal solution, and there is little difference between the optimal solution and the worst solution, the value of a 1 · P j worst − P j best is not affected by the population distribution, and this term can still get a larger weight value of w, so as to avoid falling into a small search range in the initial iteration. As the number of iterations increases, the value of a 1 · P j worst − P j best will gradually become smaller, and its influence on the weight w will decrease. If the algorithm does not get the optimal solution at this time, the design of a 2 · x maxj i − x minj i /t can play a leading role in the weight value w, which can make the algorithm search with a larger step size. The advantage of this adaptive weight adjustment method is determined by two parts. The first part changes the population when the number of iterations is too large, and the second part changes the population when it falls into the local optimal.

The Specific Procedure of the IARO Algorithm
The flow chart of the IARO is shown in Figure 3, with the steps as follows: Step 3: Calculate the weight w according to Equation (23); update parameters L , c , R , A , H , g and b ; Step 4: If

> A
, update the current rabbit position according to Equation (4); Step 5: If 1 ≤ A , update the current rabbit position according to Equation (28); Step 6: The population is updated and the best and worst rabbits found so far are updated; Step 7: If the algorithm reaches the maximum number of iterations, output the optimal solution; otherwise, return to Step 3.

Energy Shrink of the IARO Algorithm
In the iterative process of the IARO algorithm, the transition process of the energy factor from exploration to exploitation is shown in Figure 4. Step 1: Parameter initialization. Set the population size N, the maximum iteration number T, and the parameter a 1 , a 2 ; Step 2: Population initialization. The initial solution is randomly generated, the individual fitness degree {F(X i ), i = 1, 2, · · · , N} is calculated, and the best and worst rabbits found so far are recorded; Step 3: Calculate the weight w according to Equation (23); update parameters L, c, R, A, H, g and b; Step 4: If A > 1, update the current rabbit position according to Equation (4); Step 5: If A ≤ 1, update the current rabbit position according to Equation (28); Step 6: The population is updated and the best and worst rabbits found so far are updated; Step 7: If the algorithm reaches the maximum number of iterations, output the optimal solution; otherwise, return to Step 3.

Energy Shrink of the IARO Algorithm
In the iterative process of the IARO algorithm, the transition process of the energy factor from exploration to exploitation is shown in Figure 4.

Experimental Results and Analysis
To test the performance of the IARO algorithm, 5 algorithms, including IARO, ARO, PSO, ASO (atom search optimization) [36] and WOA, are compared, and 23 benchmark functions are used in the experiment [37]. The 23 benchmark functions are shown in Tables 1-3, including 7 unimodal functions, 6 multimodal functions and 10 low-dimensional multimodal functions. The IARO algorithm has been carried out using the MATLAB 9.7 (R2019b) desktop computer running Windows 10 64-bit with an Intel(R) Core(TM) i5-9400F CPU 2.9 GHz processor and 8.00 GB RAM. The initial population size of all algorithms is set to 50, the maximum number of iterations is 300, and each optimization algorithm runs repeatedly 30 times. The performance of the algorithm is analyzed and evaluated by four indicators such as the mean value (Mean), standard deviation (Std), worst value (Worst) and best value (Best) of the optimal solution found so far.

Name
Function D Range

Experimental Results and Analysis
To test the performance of the IARO algorithm, 5 algorithms, including IARO, ARO, PSO, ASO (atom search optimization) [36] and WOA, are compared, and 23 benchmark functions are used in the experiment [37]. The 23 benchmark functions are shown in Tables 1-3 including 7 unimodal functions, 6 multimodal functions and 10 low-dimensional multimodal functions. The IARO algorithm has been carried out using the MATLAB 9.7 (R2019b) desktop computer running Windows 10 64-bit with an Intel(R) Core(TM) i5-9400F CPU 2.9 GHz processor and 8.00 GB RAM. The initial population size of all algorithms is set to 50, the maximum number of iterations is 300, and each optimization algorithm runs repeatedly 30 times. The performance of the algorithm is analyzed and evaluated by four indicators such as the mean value (Mean), standard deviation (Std), worst value (Worst) and best value (Best) of the optimal solution found so far.

Name
Function Step Six Hump Camel Goldstein-Price 86

Exploitation Assessment
Since the unimodal function (F1-F7) has only one optimal value, it can be used to test the exploitation of the algorithm. The results of each indicator of the six algorithms in different unimodal test functions are shown in Table 4. The best value of each indicator has been marked in bold. Figure 5 shows the convergence curves of the five algorithms in each unimodal test function.

F6
Mean  From Table 4, for functions F1, F2, F3, F4 and F7, the 'Mean', 'Std', 'Worst' and 'Best' of the IARO algorithm are the best among the five algorithms. For function F6, the four indicators of the IARO and ARO, WOA and PSO algorithms are the best. For function F5, the 'Std' provided by the WOA algorithm is the best. The 'Mean', 'Worst' and 'Best' provided by the PSO algorithm are the best. Therefore, for the unimodal functions, the IARO algorithm is superior to other algorithms when dealing with the unimodal function.
From Figure 5, for functions F1, F2, F3, F4 and F7, the IARO algorithm has higher convergence accuracy and faster convergence speed. For function F6, the convergence ac- From Table 4, for functions F1, F2, F3, F4 and F7, the 'Mean', 'Std', 'Worst' and 'Best' of the IARO algorithm are the best among the five algorithms. For function F6, the four indicators of the IARO and ARO, WOA and PSO algorithms are the best. For function F5, the 'Std' provided by the WOA algorithm is the best. The 'Mean', 'Worst' and 'Best' provided by the PSO algorithm are the best. Therefore, for the unimodal functions, the IARO algorithm is superior to other algorithms when dealing with the unimodal function.
From Figure 5, for functions F1, F2, F3, F4 and F7, the IARO algorithm has higher convergence accuracy and faster convergence speed. For function F6, the convergence accuracy of the IARO and ARO, WOA, and PSO algorithms is the best, but the convergence speed of the IARO algorithm is slower than the other three algorithms. For function F5, the convergence accuracy and convergence speed of the IARO algorithm are lower than that of the PSO algorithm but higher than other algorithms. Therefore, for functions F1, F2, F3, F4, F6 and F7, the convergence accuracy of the IARO algorithm is optimal. For function F5, the convergence accuracy of the IARO algorithm is only lower than that of the PSO algorithm. For functions F1, F2, F3, F4 and F7, the convergence speed of the IARO algorithm is the fastest. For functions F5 and F6, the convergence speed of the IARO algorithm is in the middle of the five algorithms. From the above analysis, the IARO algorithm performs the best in Figure 5.

Exploration Assessment
Since the multimodal function (F8-F23) has multiple extreme values, it can be used to test the exploration ability of the algorithm. The results for each indicator of the five algorithms on different multimodal functions are shown in Tables 5 and 6. Figures 6 and 7 show the convergence curves of the five algorithms in each multimodal function.  From Table 5, for function F8, the 'Std' and 'Worst' of the IARO algorithm are the best among the five algorithms. The 'Mean' and 'Best' of the ARO and WOA, PSO algorithms are the best among the five algorithms. For function F9, the four indicators of the IARO and ARO, WOA, and PSO algorithms are the best. For function F10, the four indicators of the IARO and ARO algorithms are the best. For function F11, the four indicators of the IARO and ARO, PSO algorithms are the best. For function F12, the 'Mean' and 'Best' of the ARO algorithm are the best. Therefore, for the multimodal functions, the IARO shows better exploration when tackling multimodal functions.     From Figure 6, for function F5, the IARO algorithm has a faster convergence speed. For function F10, the IARO algorithm has higher convergence accuracy and faster convergence speed. For functions F8, F12 and F13, the convergence accuracy and convergence speed of the IARO algorithm are at the middle level among the five algorithms.
From Table 6, for functions F8 and F15, the four indicators of the IARO and ARO algorithms are the best among the five algorithms. For functions F16 and F17, the 'Mean', 'Worst' and 'Best' of the five algorithms are the best. The 'Std' of the IARO algorithm is the best. For function F18, the 'Mean', 'Worst' and 'Best' of the IARO and ARO, PSO, and ASO algorithms are the best. The 'Std' of the ARO algorithm is the best. For function F19, the 'Mean', 'Worst' and 'Best' of the IARO and ARO, PSO, and ASO algorithms are the best.
For function F22, the four indicators of the IARO algorithm are the best. For function, F23, the 'Mean', 'Worst' and 'Best' of the IARO and ARO, PSO, and ASO algorithms are the best. From Figure 7, for functions F14, F15, F16, F17, F18, F19 and F20, the IARO algorithm has higher convergence accuracy and faster convergence speed. For functions F21, F22 and F23, the IARO algorithm has higher convergence accuracy, and the convergence speed is significantly higher than that of the ASO and WOA algorithms.
By using 23 benchmark functions to test the performance of five algorithms, the minimum values and convergence curves of four indicators of each algorithm in different benchmark functions are obtained. Through comparative analysis, it can be concluded that the IARO algorithm has the strongest ability to find the optimal solution.

Acquisition of Experimental Data
The Zhanghewan Pumped Storage Power Station is located in Shijiazhuang City, Hebei Province, with a total installed capacity of 1000 MW. Four single-stage mixed-flow reversible pump turbine units with a single capacity of 250 MW are installed. In the power station, the rated head of the water pump turbine unit is 305 m, the rated speed is 333.3 r/min, the annual power generation is 1.675 billion kW·h, the annual pumping power is 2.204 billion kW·h, the annual power generation utilization time is 1675 h, the annual pumping utilization time is 2204 h, and the comprehensive efficiency coefficient of the power station is 0.76. After the completion of the power station, it is connected to the power grid in southern Hebei and undertakes the tasks of peak shaving, valley filling, frequency modulation, phase modulation and emergency backup in the system [38,39].
The experimental data comes from the 1# Unit data of the Zhanghewan Pumped Storage Power Station. The selected three vibration states are the stability test of power generation condition, the Karman vortex street phenomenon test and the simulated fault shutdown test. For each vibration state, three parts are selected as the sensor measuring points, and the X direction and the Y direction are arranged in coordination. The six measuring points are, respectively, the X and Y directions of the upper guide bearing, the X and Y directions of the water guide bearing and the X and Y directions of the lower guide bearing.
The experimental data of three vibration states are classified as shown in Table 7. Each vibration state is classified, that is, the power generation condition stability test No. 1, Karman vortex street phenomenon test No. 2 and simulated fault shutdown test No. 3.

Vibration Signal Processing and Extraction
The different vibration states of hydraulic units are studied experimentally, and the vibration signals are denoised and extracted by the VMD decomposition. Firstly, the center frequency method [40] is used to determine the number of decomposition layers K. Secondly, the residual index method [41] is used to determine the update step tau. In each vibration state, the number of decomposition layers K and the update step tau of the six measuring points are shown in Table 8. Finally, according to the data in Table 8, for each vibration state, the vibration signals of six measuring points are decomposed by VMD, respectively, and the corresponding IMF components can be obtained. In order to unify the number of IMF components, according to Table 8, K = 10 is selected as the number of decomposition modes extracted from each group of vibration signals; that is, IMF1-IMF10 is selected as the decomposition mode X = [X 1 , X 2 , X 3 , X 4 , X 5 , X 6 , X 7 , X 8 , X 9 , X 10 ].

Feature Vectors Selection
For each vibration state, the signals of six measuring points are averaged. The obtained mean value is substituted into the time-domain feature model. After calculation and analysis, nine time-domain features are selected as the feature vectors ] of each vibration state, which are the maximum value, minimum value, mean value, peak-to-peak value, rectification mean value, variance, standard deviation, root mean square and square root amplitude. The comparison of the nine time-domain features of the three vibration states is shown in Figure 8.
of each vibration state, which are the maximum value, minimum value, mean value, peak-to-peak value, rectification mean value, variance, standard deviation, root mean square and square root amplitude. The comparison of the nine time-domain features of the three vibration states is shown in Figure 8.

SVM Classification and Identification Model Construction
The identification of vibration states of hydraulic units is a multi-classification identification problem. The SVM [42,43] is a classification algorithm to solve binary classification problems, and it cannot be directly used to solve multi-classification problems. In this paper, the "one-to-one" (OVO) strategy is combined with SVM to form a multi-classifier. Therefore, the SVM multi-classifier design is applied to the identification of vibration signals in hydraulic units. After conducting VMD decomposition, the calculated time-domain feature vector is used as the input of the SVM multi-classifier. Combined with the IARO algorithm, the penalty parameter c and kernel function parameter g of the SVM multi-classifier are optimized. The IARO-SVM identification model is constructed, as shown in Figure 9, to identify the vibration states of the hydraulic units.

SVM Classification and Identification Model Construction
The identification of vibration states of hydraulic units is a multi-classification identification problem. The SVM [42,43] is a classification algorithm to solve binary classification problems, and it cannot be directly used to solve multi-classification problems. In this paper, the "one-to-one" (OVO) strategy is combined with SVM to form a multi-classifier. Therefore, the SVM multi-classifier design is applied to the identification of vibration signals in hydraulic units. After conducting VMD decomposition, the calculated time-domain feature vector is used as the input of the SVM multi-classifier. Combined with the IARO algorithm, the penalty parameter c and kernel function parameter g of the SVM multi-classifier are optimized. The IARO-SVM identification model is constructed, as shown in Figure 9, to identify the vibration states of the hydraulic units.

Validation of the Proposed IARO-SVM Model
The parameters of the IARO-SVM model are initialized: the initial population num ber is 30, and the maximum number of iterations is 500. The time-domain feature vector of each vibration state of the hydraulic unit are divided into the training set and the tes set, in which the first 70% of data of each group of feature vectors are used as the trainin set, and the remaining 30% data are used as the test set. All data are normalized to form sample set. After training the IARO-SVM model using the training samples, the test sam ples are input into the model. Figure 10 illustrates the identification results of the IARO-SVM model. The identifi cation accuracy of the Karman vortex and the fault shutdown signals achieve 100%, an the identification accuracy of the power generation condition achieves 93.33%.

Validation of the Proposed IARO-SVM Model
The parameters of the IARO-SVM model are initialized: the initial population number is 30, and the maximum number of iterations is 500. The time-domain feature vectors of each vibration state of the hydraulic unit are divided into the training set and the test set, in which the first 70% of data of each group of feature vectors are used as the training set, and the remaining 30% data are used as the test set. All data are normalized to form a sample set. After training the IARO-SVM model using the training samples, the test samples are input into the model.  Figure 11 illustrates the identification results of the ARO-SVM model. The identification accuracy of the Karman vortex and the fault shutdown signals achieve 100%, and the identification accuracy of the power generation condition achieves 83.33%.                 Figure 14 illustrates the identification results of the WOA-SVM model. The identification accuracy of the Karman vortex and the fault shutdown signals achieve 100%, and the identification accuracy of the power generation condition achieves 80.0%. Comparing Figures 10-14, it can be found that the identification results about the power generation condition are obviously different and the IARO-SVM model achieves the best accuracy, while other models achieve poorer accuracies. This is due to that signal differentiation corresponding to the power generation condition is the lowest, reducing the identification accuracy of the models. However, the IARO-SVM model still keeps a higher identification accuracy, indicating the LARO algorithm has a better optimization performance for the SVM.  Figure 14 illustrates the identification results of the WOA-SVM model. The identification accuracy of the Karman vortex and the fault shutdown signals achieve 100%, and the identification accuracy of the power generation condition achieves 80.0%. Comparing Figures 10-14, it can be found that the identification results about the power generation condition are obviously different and the IARO-SVM model achieves the best accuracy, while other models achieve poorer accuracies. This is due to that signal differentiation corresponding to the power generation condition is the lowest, reducing the identification accuracy of the models. However, the IARO-SVM model still keeps a higher identification accuracy, indicating the LARO algorithm has a better optimization performance for the SVM.     Figure 14 illustrates the identification results of the WOA-SVM model. The identification accuracy of the Karman vortex and the fault shutdown signals achieve 100%, and the identification accuracy of the power generation condition achieves 80.0%. Comparing Figures 10-14, it can be found that the identification results about the power generation condition are obviously different and the IARO-SVM model achieves the best accuracy, while other models achieve poorer accuracies. This is due to that signal differentiation corresponding to the power generation condition is the lowest, reducing the identification accuracy of the models. However, the IARO-SVM model still keeps a higher identification accuracy, indicating the LARO algorithm has a better optimization performance for the SVM.   To further test the efficiency of the IARO-SVM model, the 2# Unit data is selected, and the same three vibration states need to be identified: the stability test of power generation condition, the Karman vortex street phenomenon test and the simulated fault shutdown test. Except for the different units selected, the methods of data collection and classification, processing and extraction, selection of feature vectors, training and identification are the same as those of the 1# Unit data.
The parameters of the IARO-SVM model, ARO-SVM model, ASO-SVM model, PSO-SVM model and WOA-SVM model are initialized in the same way: the initial population number is 30, and the maximum number of iterations is 500. The time-domain feature vectors of each vibration state are divided into the training set and test set, in which the first 70% of data of each group of feature vectors is used as a training set, and the remaining 30% of data is used as a test set. All data are normalized to form a sample set. After training and prediction, the average identification accuracy of the three vibration states of the 2# Unit is shown in Figure 16. Inspecting this figure, it is seen that, compared with other models, the IARO-SVM model obtains the best results with an average accuracy at 96.67%, followed by the ARO-SVM model (95.56%) and the WOA-SVM model (90.00%). To further test the efficiency of the IARO-SVM model, the 2# Unit data is selected, and the same three vibration states need to be identified: the stability test of power generation condition, the Karman vortex street phenomenon test and the simulated fault shutdown test. Except for the different units selected, the methods of data collection and classification, processing and extraction, selection of feature vectors, training and identification are the same as those of the 1# Unit data.
The parameters of the IARO-SVM model, ARO-SVM model, ASO-SVM model, PSO-SVM model and WOA-SVM model are initialized in the same way: the initial population number is 30, and the maximum number of iterations is 500. The time-domain feature vectors of each vibration state are divided into the training set and test set, in which the first 70% of data of each group of feature vectors is used as a training set, and the remaining 30% of data is used as a test set. All data are normalized to form a sample set. After training and prediction, the average identification accuracy of the three vibration states of the 2# Unit is shown in Figure 16. Inspecting this figure, it is seen that, compared with other models, the IARO-SVM model obtains the best results with an average accuracy at 96.67%, followed by the ARO-SVM model (95.56%) and the WOA-SVM model (90.00%).

Conclusions
This paper proposes an identification method based on ARO with the adaptive weight adjustment strategy to SVM (IARO-SVM). By introducing the adaptive weight adjustment strategy, the algorithm can adaptively change the weight size according to the distribution of the current rabbit population. Using 23 benchmark functions, IARO is compared with ARO, PSO, ASO and WOA, and the performance of IARO is better than that of other algorithms by comparing the minimum values of the four indicators and convergence curves.
This method optimizes the parameters of the SVM multi-classifier through the IARO algorithm. Meanwhile, the VMD method is used to decompose the vibration status signals

Conclusions
This paper proposes an identification method based on ARO with the adaptive weight adjustment strategy to SVM (IARO-SVM). By introducing the adaptive weight adjustment strategy, the algorithm can adaptively change the weight size according to the distribution of the current rabbit population. Using 23 benchmark functions, IARO is compared with ARO, PSO, ASO and WOA, and the performance of IARO is better than that of other algorithms by comparing the minimum values of the four indicators and convergence curves.
This method optimizes the parameters of the SVM multi-classifier through the IARO algorithm. Meanwhile, the VMD method is used to decompose the vibration status signals of hydraulic units, and the time domain characteristic indicators of IMF components are extracted as multi-dimensional time-domain feature vectors and input into the IARO-SVM model, ARO-SVM model, ASO-SVM model, PSO-SVM model and WOA-SVM model for state identification of vibration signals. The experimental results show that the identification rate of the IARO-SVM model is 97.78%, the identification rate of the ARO-SVM model is 94.44%, the identification rate of the ASO-SVM model is 90.0%, the identification rate of PSO-SVM model is 86.67%, and the identification rate of WOA-SVM model is 93.33%. Through the comparative analysis of five models, the IARO-SVM model has obvious advantages over its competitors in the state identification of vibration signals of hydraulic units, which is 3.34% higher than the closest ARO-SVM model. For the second experiment, the IARO-SVM also achieve competitive results. Therefore, IARO has obvious advantages over other algorithms in optimizing the SVM for vibration status identification of hydraulic units.
Although the IARO-SVM model has excellent performance, it also has some limitations.
(1) The amount of experimental data is limited, so in the experimental process, sometimes there will be large deviations, which will affect the accuracy of vibration state signal classification and identification. (2) Due to the experimental conditions, the data used in the paper is only the data of pumped storage hydraulic units, and the data including other types of hydraulic units cannot be verified, the state identification for other types needs further verification. However, these shortcomings are also directions for future research. On the one hand, the amount of data can be increased and the sample set capacity can be expanded when collecting data in the future. On the other hand, the vibration state signals of other types of hydraulic units are collected to verify whether the IARO-SVM model still has obvious advantages in vibration state classification and identification for other types.